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Abstract 

Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in 
an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that 
these unstable modes have a higher growth rate than the closely related electron-scale Kelvin- 
Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the 
analytic results and further reveal the emergence of mushroom-like electron density structures in 
the nonlinear phase of the instability, similar to those observed in the Rayleigh Taylor instability 
despite the great disparity in scales and different underlying physics. Macroscopic ejoupe) fields 
are shown to be generated by these microscopic shear instabilities, which are relevant for particle 
acceleration, radiation emission and to seed MHD processes at long time-scales. 
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A fundamental question in plasma physics concerns the stability of a given plasma config¬ 
uration. Unstable plasma configurations are ubiquitous and constitute important dissipation 
sites via the operation of plasma instabilities, which typically convert plasma kinetic energy 
into thermal and electric/magnetic field energy. Plasma instabilities can occur at micro¬ 
scopic (particle kinetic) and macroscopic (magnetohydrodynamic, MHD) scales, and are 
generally studied separately using simplified frameworks that focus on a particular scale 
and neglect the other. This approach conceals the role that microscopic processes may have 
on the macroscopic plasma dynamics, which in many scenarios cannot be disregarded. It is 
now recognized, for instance, that collisionless plasma instabilities operating on the electron 
scale in unmagnetized plasmas, such as the Weibel [1] and streaming instabilities 0, play a 
crucial role in the formation of (macroscopic) collisionless shocks in astrophysical [3]-[7j and 
laboratory conditions mu- These microscopic instabilities result from the bulk interpene¬ 
tration between plasmas and are believed to be intimately connected to important questions 
like particle acceleration and radiation emission in astrophysical scenarios [iniiii]. 

Plasma shear flow conflgurations can host both microscopic and macroscopic instabilities 
simultaneously, although the former have been largely overlooked. Shear flow settings have 
been traditionally studied using the MHD framework [I2HI1!, where the Kelvin-Helmholtz 
instability (KHI) is the only instability known to develop [IH]. Only very recently have 
collisionless unmagnetized plasma shear flows been addressed experimentally [K] and using 
particle-in-cell (PIC) simulations, revealing a rich variety of electron-scale processes, such 
as the electron-scale KHI (ESKHI), dc magnetic field generation, and unstable transverse 
dynamics dMi. The generated fields and modified particle distributions due to these 
microscopic processes can strongly impact succeeding macroscopic dynamics of the shear 
flow. However, whilst both the ESKHI and the dc magnetic field formation mechanism 
have been treated theoretically, an analytical description of the transverse shear instability, 
observed in kinetic simulations, is still lacking. An established analytical description is 
crucial in order to assess its relevance in different physical scenarios. 

This Letter focuses on the shear surface instability that occurs in the plane perpendicular 
to that of the ESKHI. These new unstable modes explain the transverse dynamics and 
structures observed in PIC simulations in [miiiiEniES]. We label this effect the mushroom 
instability (MI) due to the mushroom-like structures that emerge in the electron density. 

We analyze the stability of electromagnetic perturbations in the transverse xy plane of 
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a collisionless plasma shear flow with velocity profile ho = vo{x)ez. We assume a cold 
(vth ^ Vq, where Vfh is the thermal velocity) unmagnetized plasma, and impose charge and 
current neutrality, Ueo = Tiio = Hq, v^q = ViQ (subscripts e and i refer to electron and ion 
quantities, respectively) to guarantee initial equilibrium. We use a two-fluid model where 
the relativistic equations of motion of the electron and ion fluid are coupled to Maxwell’s 
equations. We assume linear perturbations in the fluid quantities of the form f = fo + Sf 
with 5f = 5f {u E C and k E M. are frequency and wavenumber, respectively) 

with all zeroth order quantities being zero except for no(x) and Vq{x). We find, for the 
perturbed current densities, 5jx = —e(l + me/mi)no5vex,Sjy = —e(l + me/mi)noSvey,Sjz = 
—e(l -|- me/mi){nQ5vez + Snvo). Substituting in Maxwell’s equations, we can simplify the 
original set of ten coupled equations in a reduced 2x2 system: e.S'E = 0, where the 
coefficient of the tensor eij are operators of the form eij = ky,vo{x),no{x))dxm 

and hE = (SEy,SEz). In order to obtain analytical results we use a step velocity shear and 
density profile of the form Vq(x) = v~ + {v~^ — v~)'H(x) and no(x) = n~ + {n'^ — n~)'H(x) 
where H is the Heaviside function. Integrating e.hE = 0 for a: 7 ^ 0 and using the continuity 


of SEy and SE^ across the shear interface, we find solutions corresponding to evanescent 
waves: SEy^^ix) = S~Ey^z{0)e~’‘i\^^ where = c^k‘^ -\- Wpe± 7 ± - = 

e^ri^{1 rrie/rrii)/eoiTie and 7 ± = 1/a/ 1 — (n=*=/c)2. By evaluating the derivative jump of the 
electric fields across the shear interface we arrive at I.ciEo = 0 , where (iEo = (hEy(O), ^^^(O)) 
and lij = aljkj_ + an = - n;^g/ 7 o)Dj\ ai 2 = -(fc^)o/^^)(<^pe/7o)-DI^ « 2 i = ^12 

and a 22 = —1 — {k'^c^/uJ^ — /c^)D^. The dispersion relation is finally obtained 

by solving det(I) = 0. In the special case, = n“ = ho and = —v~ = ho, the growth 


rate reads 



with r = Im(a;) and Du = I /70 + k^(? 

The fastest growing mode of this unstable branch (cl^T = 0) is found at fc ^ 00 (it will 
be shown later that finite thermal effects and/or smooth velocity shear profiles introduce a 
cutoff an finite k). This limit corresponds to the maximum growth rate of the instability 
Tmax for a given shear flow Lorentz factor, and it is given by Tmax/f^pe = hb/cy^. In the 
limit, n~ S> 7 + S> 1 and v~ = 0 , relevant for astrophysical jet/interstellar medium 

shear interaction, the maximum growth rate yields Fmax/^pe-i- — l/-\/ 27 +. 
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The MI and ESKHI [I2j growth rates are compared for different shear Lorentz factors 
and velocities in Figure It is clear that the ESKHI has higher growth rates than the MI 
for subrelativistic settings. However, the MI growth rate decays with % ' , slower than the 
ESKHI, which decays with as shown in Figure 0b. Therefore, given that the noise 

sources for both instabilities are similar, the MI is the dominant electron-scale instability in 
relativistic shear scenarios. 

To verify the theoretical model and better understand the underlying feedback cycle of the 
MI, we first analyze the evolution of a single unstable mode in an electron-proton {e~p'^) 
shear flow using the PIC code OSIRIS [23l [2lj. We simulate a domain with dimensions 
20 X 5 {c/ojpeY, resolved with 40 cells per c/upe, and use 36 particles per cell per species. The 
shear flow initial condition is set by the velocity field ho = -|-0.5c for < x < and 

Uo = — 0.5c for X < La;/4 U x > 3Z/a;/4, where Lj is the size of the simulation box in the fth 
direction. The system is initially charge and current neutral. Periodic boundary conditions 
are imposed in every direction. In order to ensure the growth of a single mode, both e~ 
temperatures are set to zero, and an initial harmonic perturbation = 5vxCos{kseedy) in 
the velocity field of the electrons, with Sv^ = 10“^c and kgeed = is introduced to seed 

the mode kgeed of the instability. 

The evolution of the electron density, the out-of-plane current density Jz, and the in¬ 
plane magnetic field are presented in Figure which zooms in on the shear interface at 
X = Tx/4. The initial velocity perturbation bv^ transports electrons across the velocity 
shear gradient, producing a current imbalance in bJz- This current induces an in-plane 
magnetic field (namely bBy) that, in turn, enhances the velocity perturbation bv^ via the 
Vq X bBy force. The enhanced velocity perturbation then leads to further electron transport 
across the velocity shear gradient in a feedback loop process, which underlies the growth of 
the instability in the linear stage. The surface wave character of the fields, predicted by the 
linear theory, is also observed in Figure |^bl. 

The MI eventually enters the nonlinear phase when the growing magnetic fields become 
strong enough to significantly displace the electrons and distort the shear interface. The 
nonlinear distortion of the shear interface in the electron fluid (Figure |^a2) leads to the for¬ 
mation of electron surface current filaments (Figure [^b2). These surface current filaments 
effectively translate in a strong dc (non-zero average along the y direction) out-of-plane cur¬ 
rent structure on either side of the shear interface. Note that the protons in the background 
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remain unperturbed on these time-scales due to their inertia, and also contribute to the 
formation of the dc current. The dc current structure induces a strong dc magnetic field in 
By, as seen by the uniform field fines along the shear interface in Figures [^b2-3. The dc 
magnetic field, in turn, continues to drive the shear boundary distortion via the uq x 5By 
force, effectively mixing the electrons across the shear, and further enhancing the dc current 
structure in an unstable loop. The development of the dc magnetic field has been previously 
shown to be associated to other electron-scale shear flow processes, like the nonlinear devel¬ 
opment of the cold ESKHI and electron thermal expansion effects [IH]. Here we show that 
the MI is an additional mechanism capable of driving the dc magnetic field. The nonlinear 
distortion of the shear boundary in the electron fluid, driven by the dc field physics, ulti¬ 
mately gives rise to the formation of the mushroom-like electron density structures shown 
in Figure [^a3. Interestingly, these electron density structures are very similar to those pro¬ 
duced by the Rayleigh-Taylor instability I2SH2Z! despite the great disparity between scales 
(electron-kinetic versus hydro/MHD scales) and different underlying physics. 

The evolution of the ratio of total magnetic field energy to initial particle kinetic energy 
(es/cp), for the single mode /cgeed = 27r/5 Up^c simulation, is presented in the inset of Figure 
1^ The exponential growth associated with the linear development of the instability is 
observed for 10 < tujpe < 30, matching the theoretical growth rate. The instability saturates 
at tiOpe — 40, approximately when the size of the mushroom-like density structure is on 
the order of 27r/A:seed (Figure [^a3). Various simulations with different kseed values were 
performed, verifying the dispersion relation of Eq. (Eigurej^a). 

The inclusion of thermal effects should introduce an instability cutoff at a finite wavenum¬ 
ber ^cutoff ~ 27r/AD (Ad is the Debye length), since the thermal pressure force Fp ~ k{kBT±) 
{ks and T± are the Boltzmann constant and electron temperature perpendicular to uq), can 
overcome the electric/magnetic forces that drive the unstable arrangement of the electric 
currents at high k [2H] • Yet such effects must be incorporated through kinetic theory, particu¬ 
larly when Vth ~ Vq, since the thermal expansion of electrons may mitigate the velocity shear 
gradient, impacting the development of the instability. Interestingly, however, the thermal 
expansion is greatly reduced when considering ultra-relativistic (70 3 > 1 ) shear flows, inci¬ 
dentally the regime where the MI is most significant, allowing the use (to some extent) of the 
cold fluid treatment outlined above. When considering a plasma of temperature Tr (defined 
in its rest frame) and drifting relativistically in the 2 : direction with mean velocity Vq = /^qC, 
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most of the particles have a large Lorentz factor and thus (/3^) = J — 1, where 

f(p) is the Juettner distribution [221 ED]- The velocity dispersion in the direction of the drift 
is {(13z - l3oY) = (/5f) - /5 q ~ 0. Since + 01), we find that (130) ~ 1 - /5q = 1 / 7 q 

(exact for ^ = ksTn/meC? ^ 1), restricting the thermal expansion at high 70. The average 
particle Lorentz factor increases with Tr as (7) 7o(l + (1 + h(0)0 [3Q], where ti(0) = 3 
(3/2) for ^ ^ 1 (^ <C 1); this factor modifies the growth rate of the MI to rinax/<^pe 1 /n/ 7> 
(for the symmetric shear scenario) due to the enhanced average relativistic particle mass. 
Hence, thermal expansion effects remain negligible as long as c0 (/3^)/T^ax c/c(jpe\/( 7 ) 
[IH], which implies 70 ^ 1 . 

We have verified the role of thermal effects on the MI development in an e“e''' shear 
configuration consisting of a hot relativistic jet, with 7+ = 50 and = 0.1,1,5, and a 
cold stationary background, with 7 _ = 1 and = 10“^; we consider n~ = = Uq. We 

simulate these configurations in a domain of 200 x 100 (c/u)p00 resolved with a 4000 x 2000 
grid. As expected, the MI growth rate decreases with (Figure|^b), and is consistent with 
the temperature-enhanced average relativistic particle mass predicted analytically. The evo¬ 
lution of the electron density around one of the shear interfaces of the system is presented 
in Figure]^ for the case = 1 . At early times, k < kcutoS structures are observed in the 
electron density in Figure |^al. These small-scale structures, which are essentially surface 
current filaments, magnetically interact with each other and merge (Figure |^a2-3), forming 
larger-scale structures in a similar manner to the current filament merging dynamics asso¬ 
ciated with the current filamentation instability [3T]. The evolution of the self-generated 
magnetic field structure is also illustrated in Figure In addition, saturation of the self¬ 
generated magnetic field (Bgat) is achieved when F^ax ^ ojb = based on 

magnetic trapping considerations [H 2 ] . 

We have also considered the effects of a smooth shear in the development of the MI. We 
have performed 2 D PIC simulations with velocity profiles of the form vo{x) = vq tanh(x/Ly), 
where Lj, is the shear gradient length, which reveal the persistent onset of the MI at gradient 
lengths up to Ly = c/ujpi ^ c/cUpe (where uipi is the ion plasma frequency). The MI can thus 
be of relevance to the physics on the ion length/time scale. It is found that the introduction 
of the finite shear length Ly also leads to a fastest growing mode at a finite /Cmax — 27r/L„. 
The growth rate F^ax of the MI decreases with increasing L„, but remains higher than the 
ESKHI for the same Ly in relativistic scenarios. 
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Interestingly, due to its electromagnetic nature, the MI is found to operate in the absence 
of contact between flows, i.e., in the case of a finite gap Lg separating the shearing flows 
(Figurej^, highlighting the different nature of the MI compared to the (bulk) two-stream and 
current filamentation instabilities. This setting is closely connected to the work explored in 
[33] on the development of optical instabilities in (subrelativistic) nanoplasmonic scenarios, 
consisting of shearing metallic slabs separated by a nanometer-scale gap; the development 
of such instabilities results in an effective non-contact friction force between slabs [311 EH] • 
The role of the MI, however, was overlooked since only subrelativistic configurations were 
considered; the transverse MI mode will be predominant in the relativistic regime. The new 
MI modes are found by taking into account the new boundary conditions imposed by the 
gap. The new eigenmodes of the system have a more complex spatial structure; the surface 
waves peak at the flow boundaries, evanesce with k\_ |no=neo k±_ |„(,=o in the plasma 
and vacuum regions, respectively, and couple in the gap. Therefore, the interaction between 
flows across the gap is strong as long as |no=o"^ 1- Simulations show that the most 

unstable mode is found at ~ 27r/Lp and that the instability growth rate decreases as 
the gap is broadened, T ^^p , ~ Tmax exp(—L^cUpe/cv^). 

In conclusion, we have described a fast-growing electron-scale instability (MI) that devel¬ 
ops in unmagnetized shear flows, and that explains the structures observed in HZmsiEolEg. 
We have shown via analytic theory that the MI grows faster than the closely related ESKHI 
(that operates in the perpendicular plane) in relativistic shears. These analytical results 
support the findings reported in [El, where 3D PIC simulations of cold, unmagnetized 
shear flows were investigated, showing the simultaneous operation of the ESKHI and MI, 
and where it was found that the MI dominates over the ESKHI in relativistic scenarios. 
Eurthermore, we have presented 2D PIC simulations that give insight into the formation of 
the mushroom-like electron density structures in the nonlinear stage of the MI. Finally, the 
MI was shown to be robust against a variety of non ideal (realistic) shear conditions: (i) 
operating under relativistic thermal effects, (ii) different densities between shearing flows, 
(iii) smooth velocity shear profiles [L^ c/ujp^), and, quite surprisingly, (iv) even in the ab¬ 
sence of contact between flows. Relativistic shear flow conditions may be reproduced in the 
laboratory by propagating a globally neutral relativistic e“e+ beam [36] in a hollow plasma 
channel [371 Ei, allowing experimental access to the MI and shear flow dynamics on the 
electron scale. This will be explored and presented elsewhere. In the astrophysical context. 
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the unexpected ability of these microscopic shear instabilities (MI and ESKHI) to generate 
strong (equipartition) macroscopic in the case of e~ shears) fields from initially 

unmagnetised conditions can directly impact particle acceleration and radiation emission 
processes, and can seed the operation of macroscopic (MHD) processes at later times (e.g. 
magnetic dynamo). 
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FIG. 1. Maximum growth rates of the MI and ESKHI versus (a) flow velocity Fq and (b) flow 
Lorentz factor 70 , highlighting the dominance of the MI (ESKHI) in relativistic (subrelativistic) 
flows. The red points represent results of 2 D PIC simulations. 
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FIG. 2. Evolution of a single unstable mode (A:seed = 2n/5LOpe/c) of the MI in an e“p+ shear. 
Frames (a) and (b) show the electron density and the out-of-plane current density J^, at times 
(1) t = 16/wpe, (2) t = 26ju!pe and (3) t = 38/ojpe. The evolution of the self-consistent in-plane 


magnetic field is also displayed in Frames (b). 
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FIG. 3. a) MI growth rate versus wavenumber for an e~p+ shear flow with = —v~ = Fq = 0.5c 
(solid line); red dots represent results of 2 D PIC simulations. The inset shows the evolution of 
^bI^p for tho single mode simulation with /cgeed — ^pejc- b) Finite temperature effects on 

the growth rate of the MI triggered by a hot relativistic e“e+ jet with 7 + = 50, shearing with a 
cold stationary plasma with 7 _ = 1 and = n~ . Dashed lines represent the theoretical slope 
rmax/‘^pe+ = 17 ^ 2 ( 7 +) for each case. 
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FIG. 4. Electron density evolution in relativistic e“e+ shear, with 7 + = 50 and 7 _ = 1, at (al) 
t = SO/cjpe, (a 2 ) t = IGO/cjpe and (a3) t = dOO/cj^e- Small-scale current filaments are excited at 
early times, merging into large-scale structures as the instability develops. Vector field represents 
self-generated magnetic field structure. 
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FIG. 5. MI development in the case of a finite gap between flows 
Fo = ^2/3 c. Frames a) and b) reveal the electron density at times t 
respectively. 


with Fgap = 5 c/ujpe and 
= O/cjpe si^nd t = 72blujpe^ 
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